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We introduce a new time-domain method for computing the self-force acting on a scalar 
particle in a Schwarzschild geometry. The principal feature of our method consists in the 
division of the spatial domain into several subdomains and locating the particle at the 
interface betweem two them. In this way, we avoid the need of resolving a small length 
?H , scale associated with the presence of a particle in the computational domain and, at the 

same time, we avoid numerical problems due to the low differentiability of solutions of 
equations with point-like singular behaviour. 
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1. Motivation 

One of the main sources of gravitational radiation for the future space-based 
gravitational-wave observatory LISAE are the capture, and posterior inspiral, of 
stellar-mass compact objects (SCOs) into massive black holes (MBH) located at 
galactic centers. Since the masses of interest for the SCO are around m = 1 — 10 2 M Q , 
and for the MBH are in the range M = 10 4 — 10 7 M Q , the mass-ratio for these systems 
is fi = m/M ~ 1CP 7 — 10~ 2 . For this reason, they are called Extreme-Mass-Ratio 
Inspirals (EMRIs). During the inspiral phase, an EMRI losses energy and angular 
momentum via the emission of gravitational waves (GWs). LISA will be able to 
detect GW signals of 10 — 10 3 EMRI/yr up to distances with z < lP These signals 
will be hidden in the LISA instrumental noise and in the GW foreground produced 
mainly by compact binaries in the LISA band. Thus, in order to extract the EMRI 
signals we need a very accurate theoretical knowledge of the gravitational wave- 
forms. The main difficulty in producing those waveforms is the description of the 
gravitational effects of the SCO on its own trajectory. These effects produce devia- 
tions in the motion of SCO, which is not longer a geodesic around the MBH, which 
can be pictured as the action of a local force, the self-force. Here, we review the 
results of recent worlP^ where a new time-domain technique for the computation 
of the self-force has been proposed. 
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2. Top tips for an efficient time-domain self-force computation 

Due to the complexity of the gravitational EMRI problerrpE' we have developed our 
new techniques for self-forse computations using a simplified model that contains 
all the ingredients of the gravitational case. It consists of a charged scalar particle 
(the SCO) in a circular motion around a non-rotating MBH. This provides a good 
test bed to test our techniques before extending them to the gravitational case. 

The equations of a (scalar, electromagnetic, gravitational) field on a 
Schwarzschild background inherit the spherical symmetry of the geometry. There- 
fore, we can decompose the field in harmonic modes, eliminating the angular de- 
pendence, so that each of them satisfies a decoupled 1+1 wave-type equation. In 
contrast with the behavior of the full field, each harmonic mode turns out to be 
finite at the particle location, which is very useful to regularize the field mode by 
mode using the mode-sum regularization scheme.^® Then, we can obtain a regular 
field by adding all the regularized harmonic modes, from which we can compute 
the self-force acting on the particle. It is then very important to develop efficient 
techniques to compute precisely the harmonic modes near or at the particle location 
so that we can estimate very precisely the self-force via the mode-sum scheme. 

Recently, we presented 3 ^ a new time-domain technique to solve the wave equa- 
tions for the harmonic modes in an efficient and precise way. It consists in a multi- 
domain framework where the particle is always located at the interface between two 
subdomains. This setup has two important advantages: (i) We do not need to re- 
solve a small scale associated with the presence of a point-like singularity inside the 
computational domain and (ii) to avoid the negative effects in the numerical com- 
putations of the low differentiability of the solution. The main idea is that the wave 
equations in the different domains are source-free and then do not see point-like 
singularities, which ensure good differentiability properties and hence, good numer- 
ical convergence properties. The way the different domains are communicated is 
through analytical junction conditions dictated by the field equations themselves.^ 

Regarding the numerical implementation, we perform the spatial discretiza- 
tion using a Chebyschev-Lobatto Pseudo spectral Collocation Method^ (PCM). We 
evolve a first-order system of equations, obtained from the reduction of the wave- 
type equations, that allow us to impose the junction conditions on the characteristic 
fields of this system. In practice, this is imposed via the penalty method^^ which 
drives the system dynamically to satisfy the junction conditions pEl The convergence 
properties of the PCM are very sensitive to the smoothness properties of the so- 
lution, which gives more importance to the multi-domain techniques we are using. 
In this sense, it has been showrPEl that our methods are able to resolve with pre- 
cision the field on the particle location with a reasonably low computational cost, 
showing that indeed the method is well suited for time-domain computations of the 
self-force. Up to now, we have done calculations for circular orbits. Comparing our 
results with others in the literature) 13 * 14 ^ we have found that they agree with a high 
degree of precision even when we use a relatively low number of collocation points. 



January 26, 2010 14:20 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



3 

For instance, at the last stable circular orbit (r — 6M), our computation of the 
radial derivative of the field, the only one that requires regularization in the circular 
case, coincides with the values obtained in other time-domain and frequency-domain 
calculations with a relative error of the order of 0.2%. Our calculations use between 
12 and 24 subdomains (for the discretization of the radial direction in terms of the 
tortoise coordinate) and 50 collocation points per domain. The average time for a 
full self- force calculation (which involves the calculation of 231 harmonic modes) in 
a computer with two Quad-Core Intel Xeon processors at 2.8 GHz is always in the 
range 20-30 minutest These calculations can be further optimized by distribut- 
ing the subdomains and collocation points so that the resolution is adapted to the 
physical problem. The calculations can be easily parallelized, either by spreading 
the work of the harmonic modes or that of the subdomains. 

Looking at the future, we are currently extending these techniques to eccentric 
orbits. This has required some modifications to keep the particle fixed at a node 
between subdomains, and results of the calculations will be published elsewhere.^ 
The next step will be the gravitational case, where the challenge comes from the fact 
that each harmonic mode is described by a set of coupled 1+1 wave type equations. 
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